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Abstract 

We develop a simple diagrammatic low-energy boson exchange strong-coupling 
model for an antiferromagnetically correlated polaronic metal. The model is based 
on a simplified, local anharmonic oscillator representation of the highly anharmonic 
Born-Oppenheimer lattice potential which typically arises in a finite density po- 
laron system, as a result of polaronic self-localization. This Born-Oppenheimer 
potential exhibits a large manifold of nearly degenrate local potential minima, 
and hence a new class of low-energy excitations, associated with tunneling pro- 
cesses between the multiple Born-Oppenheimer wells. In the present paper, we 
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Q\ . der to study the what happens when conduction electrons are coupled to these 
anharmonic tunneling excitations via conventional electron-phonon interactions. 
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Important low-temperature features of such a model are (i) a large polaronic mass 
enhancement Z and a strongly narrowed, low-energy, near- Fermi- level van Hove 
resonance peak in the interacting density of states at temperatures T well be- 
low the lattice tunneling energy scale Q t ', (ii) a rapid thermal suppression of Z 
and of the van Hove resonance which occurs when the temperature T becomes 
comparable to the lattice tunneling excitation energy scale Qf, (hi) a broad inco- 
herent background in the single-particle spectra \G"(k, u)\ near the Fermi energy; 
(iv) strongly flattened quasi-particle bands over wide regions of the Brillouin zone 
near the van Hove points; and (v) strong isotopic mass dependence of the elec- 
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tronic mass enhancement Z at temperatures T well below Q t which, along with 
Z, becomes rapidly suppressed for T comparable to Cl t . 

Surprisingly, in spite of the large low-T mass enhancment, the anharmonic 
lattice fluctuations in the high temperature regime T ~ f^, give rise only to a 
moderate, essentially T-independent quasi-particle damping. This damping does 
not suppress coherent quasi-particle propagation in a wide, physically relevant 
parameter regime. We thus argue that, in a "weakly bound" finite-density polaron 
system, a novel polaron scenario may be realized where only the lattice motion, but 
not the quasi-particle propagation become incoherent in the high-T regime. This 
scenario is in contrast to the conventional small-polaron strong-coupling picture, 
where the quasi-particle transport inevitably becomes incoherent at temperatures 
exceeding the polaron tunneling energy scale. 

We also discuss possible implications of the anharmonic tunneling model for 
the cuprate superconductors. The predicted isotopic mass dependence of the low- 
T (T < Q t ) electron quasi-particle mass enhancement provides a unique signature 
for the involvemnt of anharmonic tunneling fluctuations. We suggest possible low 
temperature isotope experiments to explore this question in the cuprates. 

PACS 74.20.-z., 75.10.-b., 63.20.Kr, 71.30.+h, 71.45.Lr 
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I. INTRODUCTION 

Hubbard-type strongly correlated electron models 1 can adequately describe 
the low-lying charge and anti-ferromagnetic (AF) spin excitations in the undoped, 
insulating cuprate systems. 2-5 However, in the doped systems, the low-energy 
physics may be substantially affected by coupling to the lattice degrees of freedom, 
since, in the presence of AF spin correlations, already a rather weak electron- 
phonon coupling can cause polaronic self- localization of the dopant-induced charge 
carriers. 6-10 As a consequence, large-amplitude anharmonic lattice displacement 
fluctuations occur, resulting from polaron tunneling motion. 6-7 In the cuprate 
high-T c systems, there is indeed evidence for polaron formation and for large- 
amplitude, anharmonic lattice tunneling fluctuations of possible polaronic origin. 11 
However, the observed quasi-particle mass enhancements Z ~ 2 — 4 in the normal 
state of the doped, metallic cuprates 12-19 fall far short of the magnitudes Z ~ 
10 — 10 3 (or larger) which one would naively expect to observe in a polaronic 
system. 

In the present paper, we propose a new modeling approach towards the 
normal-state electron quasi-particle dynamics of an antiferromagnetically corre- 
lated polaronic metal where the essential features of the polaronic lattice 
dynamics 6-7 and AF spin fluctuation exchange 20-29 are combined into a gener- 
alized Migdal-type diagrammatic treatment. The polaronic lattice anharmonicity 
does then indeed lead to large polaronic mass enhancement Z when the tempera- 
ture T and excitation energies u are well below a certain, generally small lattice 
tunneling energy scale Clf 

However, for T ~ fi t , the polaronic contribution to Z is suppressed and Z 
is dominated by the smaller spin fluctuation contribution. The polaronic lattice 
fluctuations contribute then only to the quasi-particle damping, in a manner quite 
similar to - and experimentally indistinguishable from - disorder scattering. Sur- 
prisingly, in spite of the large low-T mass enhancement, the high-T quasi-particle 
damping caused by the anharmonic lattice fluctuations is quite modest and the 
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electronic transport dynamics retains essentially coherently propagative character. 
This suggests the possibility of a heretofore unexplored "weakly bound" polaron 
scenario wherein polaron formation does not inevitably lead to complete incoher- 
ence at temperatures exceeding the polaron tunneling energy scale. Our treatment 
thus suggests a novel approach towards dense polaron systems which bridges the 
gap between the conventional Fermi-liquid-based Migdal theory of delocalized, har- 
monic electron-phonon systems 30 ' 31 and the strong-coupling methods developed 
originally for the strongly localized, single-electron small-polaron problem. 32-37 

The remainder of the paper is organized as follows: In Section II, we outline 
the basic underlying ideas of polaronic anharmonicity in the lattice dynamics of 
a strongly correlated electron-phonon system; we propose a simple exactly nu- 
merically solvable model for the polaronic Born-Oppenheimer potential; and we 
discuss the effective electron-electron interaction Vap which is mediated by ex- 
change of the anharmonic tunneling excitations arising from the lattice quantum 
dynamics within such an "anharmonic phonon" potential. We then combine this 
lattice-mediated electron interaction Vap with a simple electron-electron model in- 
teraction Vsf describing AF spin fluctuation exchange. The full model interaction 
V = Vap + Vsf then serves as our input into a self-consistent diagrammatic single- 
electron self-energy calculation, which we formulate at the level of a Migdal-type 
single-exchange approximation for a single-band conduction electron system on a 
two-dimensional (2D) square lattice. From the self-energy, we extract important 
physical quantities such as the electron quasi-particle dispersion Er, and and the 
quasi-particle mass enhancement Z, the single-electron spectral function G"(k,u)) 
and the the momentum-integrated density of states p(u>). 

In Section III, we discuss basic physical principles governing the polaronic 
mass enhancement in our diagrammatic model. We show, by analytical arguments, 
how, within our diagrammatic model, one can recover, at least qualitatively, the 
central results for polaronic mass enhancement which are conventionally derived 
within a non-perturbative electron-phonon strong-coupling theory. Important im- 
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plications for a potential observable isotopic mass dependence in the electron mass 
enhancement are discussed in the anharmonic tunneling exchange model and con- 
trasted with a corresponding harmonic phonon exchange model. Based on self- 
consistent numerical self-energy calculations, we then discuss the temperature (T-) 
dependence of the polaronic self-energy and its implications for the T-dependence 
of the quasi-particle mass enhancement Z and of the quasi-particle damping. We 
then consider the effects of anharmonic tunneling fluctuation exchange on the mo- 
mentum integrated density of states and, in particular on the van Hove singular- 
ity. Numerical results for the single-particle spectral functions and for the overall 
quasi-particle band dispersion in the full Brillouin zone will also be discussed. 

Lastly, in Section IV, we present a summary of our results, with emphasis on 
potential experimental implications of our model calculations for the cuprate high 
T c materials. We conclude with an outline of the generic features of our model 
which are likely to be of general physical significance and likely to be recovered in 
a more realistic polaron model. 

II. LATTICE TUNNELING FLUCTUATION EXCHANGE MODEL 

A. Polaronic anharmonicity 

As a simplest model, we start from a two-dimensional (2D) antiferromag- 
netically correlated Hubbard-type single-band electron system on a 2D square 
lattice, 1 ~ 5 ' 20-29 coupled to an Einstein optical phonon system via a Holstein 
electro n-phonon interaction 6 ~ 10 ' 32 ~ 37 



j 

Here, C denotes the deformation potential coupling constant, nj = 0, 1, 2 the elec- 
tron occupation, and uj the local Einstein oscillator displacement at lattice sites 




(1) 
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j = l...N. Physically, Uj could represent, for example, a collective coordinate for 
local breathing-type displacements of planar or apical oxygens in the cuprates. 

In the absence of EP coupling, the lattice dynamics is governed by a bare 
harmonic restoring potential Wp({uj}) corresponding to a non-interacting phonon 
Hamiltonian 

H r = Y.m p2 i +Wp{{u > }) (2) 

3 

where M is the atomic or molecular mass associated with each displacement de- 
gree of freedom Uj and pj = —ihd/duj is the conjugate momentum of Uj. In 
the simplest microscopic strongly correlated electron-phonon model, the Holstein- 
Hubbard model 6 ' 7 ' 9 ' 32-37 Wp is given by an Einstein potential 

Wp({u j }) = J21ku 2 j (3) 

3 

with a bare, harmonic restoring force constant K. To be specific, we will, for the 
following discussion, adopt the Holstein-Hubbard model as our underlying micro- 
scopic model Hamiltonian. The basic ideas are more general, however, and can, 
in principle, be extended to more complicated electron-phonon systems exhibiting 
polaronic self-localization. The EP coupling strength can then be quantified by a 
characteristic energy 

r 2 

E P = ^ (4) 

which measures the strength of the phonon-mediated on-site attraction, as well 
as the ingle-polaron binding energy in the ionic (zero-bandwidth limit) of the 
Holstein-Hubbard model. 

In the nearly | -filled Holstein-Hubbard electron system with sufficient EP 
coupling strength Ep, a dopant-induced carrier of, e. g., hole-type, rij = 0, self- 
localizes. 6-9 That is, by locally distorting the lattice by some Auj = d tl the carrier 
creates an attractive electron potential well, of depth 



A P = Cd t 
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, and thus lowers its energy by forming a bound state. 6 ~ 9 ' 32-37 The inter-site 
tunneling motion of such polaronic carriers causes the displacement Uj at each site 
j to fluctuate between two distinct equilibrium positions, u±, with \u + —u_ \ = d t . 

As explained at length in Refs. [6,7,38], this physical picture can be fully cap- 
tured within the framework of the Born-Oppenheimer approximation, that is, in 
terms of a renormalized lattice potential W({uj}). 6,7 ' 38 In the full microscopic 
treatment, based on the Holstein-Hubbard model, 6 ' 7 this Born-Oppenheimer po- 
tential is given by the total ground state energy of the Oth order adiabatic Hamil- 
tonian H a d({uj}) which depends parametrically on the fixed, c-number lattice 
displacements {uj}. H ad ({uj}) consists of the bare potential Wp({uj}), the EP 
coupling term HE({uj}) and the purely electronic part He which is independent 
of the displacements Uj. Tn the case of the Holstein-Hubbard model, He is just 
the standard Hubbard Hamiltonian. 1 ' 2-5 Detailed microscopic studies along these 
lines, based on analytical considerations and on numerical many-body total energy 
calculations on finite clusters have been carried out. 6 ' 7 ' 9 

These microscopic studies reveal that the transition or cross-over from the 
delocalized carrier regime to the self-localized carrier regime is fundamentally ac- 
companied by a qualitative change in the character of the renormalized lattice 
potential W. In the delocalized carrier regime, the low-energy lattice dynamics 
is governed essentially by one lowest, absolute potential minimum in W({uj}). 
Quadratic expansion of W({uj}) around its minimum, thus leads again to an es- 
sentially harmonic lattice dynamics where the primary renormalization effect of 
the EP coupling is to modify the harmonic restoring force constants, relative to 
Wp({uj}). 

In the self-localized carrier regime, on the other hand, W({uj}) acquires a 
highly degenerate, anharmonic, multiple-minimum low-energy structure wherein 
different local minima of W correspond to different possible real-space config- 
urations of the self- localized carriers. The polaron tunneling dynamics in this 
Born-Oppenheimer picture can then be described in terms of lattice tunneling 
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processes between "adjacent" VF-minima. The low-energy structure of the Born- 
Oppenheimer energy surface is thus highly complex, especially, at a finite polaron 
density, involving, in general highly anharmonic, highly non-local multi-site in- 
teractions betwen the Uj coordinates. The low-energy lattice tunneling dynamics 
within this degenrate multi-well structure is further complicated by the presence 
of non-trivial e l7F Berry phase factors which arise from adiabatic electron ground- 
state wavefunction overlaps during the adiabatic tunneling motion of the lattice 
between different local W-we\\s. 7 Trying to calculate, or even simulate, the low- 
energy lattice dynamics in the self- localized regime, from a "realistic", that is, 
microscopically based Born-Oppenheimer potential does not appear to be feasi- 
ble at the pesent time, especially in the physically most interesting case of finite 
polaron carrier densities. 

We are therefore proposing to simplify the polaronic lattice dynamics problem 
by replacing the microscopic lattice potential with a drastically simplified model 
potential. This simplified potential does not, by any measure, capture all the 
intricate details of the microscopic potential. Nevertheless, as we will show, it 
retains some of the essential physical features of the polaronic lattice dynamics. 
Furthermore, this model potential can serve as the starting point for more realistic 
modeling approaches which may capture some of those physical aspects, such as 
non-local multi-site anharmonicity, which our present model fails to include. 

Specifically, we propose to model the polaronic lattice potential by the follow- 
ing ansatz 



The single-site potential w(uj) is a local (on-site) double- well with minima at 



The two potential minima in w(uj) are separated by a tunneling barrier of height 
Ab and width dt. 





u± = ±d t /2 . 



(7) 
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Eq. (6) captures some of the essential physics of the polaronic lattice dynamics, 
namely, the fact that each Uj undergoes large-amplitude fluctuations between two 
distinct equilibrium positions u±. Neglected in this simple model potential are are 
all the anharmonic inter-site correlations described above within the microscopic 
framework. 6 ' 7 

The essential new feature of the polaronic lattice dynamics, which is absent in 
harmonic phonon systems, is a new class of large-amplitude, low-energy inter-well 
tunneling excitations, which exist in addition to the familiar, small-amplitude, 
intra- well phonon-like excitations. In our model, Eq. (6), these two types of ex- 
citations are represented by the excitation energies fi t and Vt' h of, respectively, 
the 1st and 2nd excited state, measured from the groundstate, in the single-site 
double- well w(u). For reasonable parameters, we find 

n t « n' h (8) 

where Vt' h is comparable to the single-well harmonic phonon energy which can be 
estimated by 

where the single- well hamonic restoring force 

K ^J*L w(u±)=40 ^ (10) 

is determined by the curvature of the potential at the single-well minimum and M 
is the atomic/molecular mass associated with the local oscillator degree of freedom 
u. In a full microscopic Holstein-Hubbard calculation, 6 is generally comparable 
to the bare (C = 0) harmonic phonon energy fi^o = h(K /M) 1 / 2 . 

B. Fluctuation exchange potentials 

By exchange of anharmonic lattice excitations an effective electron-electron 
interaction Vap is mediated. 39-43 As in standard harmonic phonon exchange 
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models, 30 ' 31 this anharmonic phonon exchange potential VAp(q, iw) can be ex- 
pressed as a product of electron-lattice coupling constants [i.e C 2 in our model, 
Eq. (1)] and appropriate lattice displacement correlation functions. 39-43 That is, 
in our model, following Refs. [39-43], 

V AP (q,iuj) = -C 2 [\ T ^Te^ T -^-^( Uj (T) A pu e (0)) A p . (11) 
J ° N iJL 

for an electron-electron momentum-energy transfer (q, iui) where iu = 2n L0 TiT is 
the even Matsubara frequency with integer n w . N denotes the number of lattice 
sites on an N = L x L 2D square lattice lattice with periodic boundary conditions, 
fj and fi are the position vectors of lattice sites j and £, respectively, (3 = 1/T is 
the inverse temperature. The notations (...)ap an d Uj(r) A p in Eq.(ll) indicate 
thermal averaging and, respectively, evolution in imaginary time r, according to 
a Born-Oppenheimer effective lattice Hamiltonian given by 

= E + w ^)) = E v • ( 12 ) 

3 3 

For our simple local model potential, Eq. (6), Hap describes N uncoupled local 
double- well oscillator degrees of freedom Uj, each subject to the local Hamiltonian 

K = t^P 2 + w(u) , (13) 

which acts only on a single local degree of freedom u. 

Because the simple on-site coupling, Eq. (1), and local model potential, 
Eq. (6), V A p = V A p(iuj) becomes ^-independent in our model, and can be ex- 
pressed in terms of the local (single-site) eigenstates of h u . 39-43 Namely, written 
in terms of its spectral representation, 

J-oo n u'-iu) 

where the spectral function on the real-c/ axis, V a ' p (uj') = \\V A p{oj' + i0 + ) — 
V A p(uj' — ^0 + )], is given by 

Vap(w') = -7T V |<^l*}| 2 [S(Vk -V,- - 5(77. - rj L + "')] • (15) 

T7? z h 



Here, ip L (u) denotes the normalized eigenstates of h u with eigenenergy rj t , labeled in 
ascending order, beginning with the groundstate t = 0. Also, zu = J2 L ex P(~PVt) is 
the corresponding single-site partition function. V^ p can thus be easily evaluated 
by solving numerically the single-site Schrodinger equation for h u , to obtain ip t 
and rj L . By carrying out the t- and k- summation in Eq. (15) numrically up to, 
say, the 5-th excited state, i, k < 5, one captures over 99% of the total spectral 
weight in V^ P , assuming physical parameter values given below. 

We augment Vap by a simple model for the electron-electron interaction po- 
tential, VsF(q,iw), mediated by the AF spin fluctuations. 20-29 The model poten- 
tial is defined for the 2D square lattice via its spectral spectral representation 

V S F{q,w)= ^7 : — 16 



oo 

with a spectral function 44 of the form 24 

vgM«>) = -g 2 sY,nf -&)*"{<*) ■ (17) 

Q* 

A simple continuous, ^-independent linear spin fluctuation spectrum is assumed, 
with 24 

a>» = ^e(n s - H) . (is) 

extending up to a typical spin fluctuation cutoff O s . The ^-dependence is assumed 
to be Lorentzian 24 

A S F 



F(q-Q*) = ~ 2 (19) 



|<f-Q*| 2 + « 2 



with a half-width kaf of order of the inverse AF correlation length. The Q* 
summation in Eq. (17) is over the four equivalent AF wavevectors 

Q* = (7T, 7T), (-7T, 7T), (-7T, 71"), (-7T, -7r) . (20) 

The overall normalization factor Ag^? is chosen so that 

— 1 f°° 

9s = ^J2j o d " V s' F (q,") (21) 



where gl is the spin fluctuation coupling constant and the ^-summation is only 
over the 1st Brillouin zone of the 2D square lattice with lattice constant a. Also, 
Eq. (17) is understood to apply only to g-points within the 1st Brillouin zone, with 
periodic continuation beyond the 1st zone boundaries. 

We should emphasize here again that the foregoing simplified model potentials 
can reproduce only the most qualitative physical features of the order of magni- 
tude, of the T-variation and of the ^-dependences of various physical quantities 
discussed below. For instance, no T-dependence is built into the spin fluctuation 
exchange potential. Likewise, the spectral and thermal properties of the polaronic 
lattice dynamics has been severely simplified here, to retain physical features such 
as large-amplitude tunneling fluctuations and a low tunneling (Q t ) energy scale, 
while neglecting, e. g., all effects from inter-site lattice couplings which, surely, will 
alter the details of the T- and w-dependence of V^ P . Therefore, these potentials 
cannot and should not be expected to reproduce, say, details of the quasi-particle 
spectral shapes or of the T-dependence of quasi-particle lifetimes. 



C. Self-energy calculation 

From the total interaction V = Vap + Vsf, we obtain the single-electron self- 
energy E in a self-consistent Migdal single-exchange approximation. 30 ' 31 ' 44 In the 
Matsubara frequency domain, this single-exchange diagram for the single-electron 
self-energy is given by 30 ' 31 ' 44 

E(jfe, iu) = -^ V (k ~ k', %v - ii/)G(k', iv) (22) 

k',iu' 

where k and kl denote electron momenta and iv = (2n u + 1)ttT and ^ / (2n' t , + 1)ttT 
are odd Matsubara frequncies with integer n v and n' v . 

The fully dressed electron Green's function in Eq. (22) is given by the Dyson 
equation 

G$M = - 1 r , - (23) 

%v — eg — iv) 
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where eg is the non-interacting single-electron conduction bandstrucure. We are 
adopting a single-band tight-binding model with 1st neighbor hybridization t only, 
again defined on the 2D square lattice appropriate to a single Cu02 layer in the 
cuprates. Hence, for k = (k x ,k y ), 

eg = — 2t[cos(ak x ) + cos(ak y )] — n (24) 

with the chemical potential \i absorbed in eg. 

The basic idea underlying this diagrammatic approach is that, at the level of 
the Born-Oppenheimer approximation, all renormalizations of the phonon system 
due to EP coupling are, in principle, contained in the effective lattice potential W. 
Hence, the displacement correlation function in Eq. (11) represents, in principle, 
the fully dressed phonon propagator and Vap must, in this context, be regarded 
as the already fully screened phonon-mediated electron-electron interaction po- 
tential, which can be directly inserted into the Migdal single-exchange diagram, 
without further renormalization. Beyond the well-controlled Born-Oppenheimer 
approximation, there are thus two primary further approximations being made 
here: firstly, the replacement of the quite complex polaronic multiple-well lattice 
potential by the simple local double- well model, Eq. (6), and, secondly, the neglect 
of higher-order vertex corrections, in the Migdal approximation. Of these two, the 
former is likely to be the more severe one; the latter can actually be roughly jus- 
tified for the model parameter values employed below, based on standard order of 
magnitude estimates for the leading order vertex correction. Another essential - 
and likely uncontrolled - approximation is of course the treatment of AF spin fluc- 
tuation exchange within the Migdal approximation. 45 This constitutes a common, 
fundamental problem in all presently existing diagrammatic AF spin fluctuation 
exchange theories. We have, in that regard, no further progress to offer in the 
present paper. 

By standard analytical continuation techniques, 44 Eq. (23) can be recast into 
the real-frequency domain to yield the self-energy spectral function E" on the 



real-w axis, expressed in terms of the spectral functions of V and G, namely 
oj) = -[£(£, oj + i0+) - E(fc, w - i0+)] 



fe' 

x [6(w - oj') + 1 - /(a/)] • (25) 



Here, /(a;) = l/(e^ w + 1) and 6(u;) = l/(e (3uJ — 1) denote the Fermi and Bose 
factor, respectively, and the /^-integral is over the 1st Brillouin zone of the 2D 
square lattice. The spectrum of the single-electron Green's function in Eq. (25) is 
given by 

G"(k, oj) = - [G(k, tv + iO+)- G(k, oj - i0+)] 

** W > (26) 



[w - eg - E'(fc, oj)] 2 + [E"(fc, oj) + 0+] 2 



and E' is obtained from the self-energy spectral function E", via the Kramers- 
Kronig relation 

Z'{k,uj) = l[E'(k,u + iO+) + J:'(k,u-iO + )] = J>— ^1^1 (27) 
2 J TV oj — OJ 

where If denotes the principal value integral. 

Eqs. (25)- (27) are solved iteratively, starting from the 0-th order ansatz for the 
self-energy Eq = 0. To evaluate E" for N — > oo numerically in the real-frequency 
domain, Eq. (25), with high ^-resolution, we use a modified, 2D adaptation of 
the so-called "tetrahedron method" 46 to carry out the fc'-integral in Eq. (25). 
Specifically, to handle the k' integration over the nearly singular quasi-particle 
peak structure of G" efficiently, the inverse of the electron Green's, from Eq. (23), 
that is 1/G(k, oj ±i0 + ) [not G{k, oj ± i0 + )] and, separately, the (^-dependent factor 
of the interaction potential, V", from Eqs. (17), (18), are linearly interpolated over 
small triangular /c-elements. Upon, expressing G" = ImG in terms of the G, so 
interpolated, the /c'-integral over each such triangular element can be calculated 
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analytically. To construct the triangular interpolation grid, the Brillouin zone is 
covered with an N = L x L grid of equal square elements where each square is then 
subdivided along a diagonal into two equal triangles. To preserve the 4-fold point 
symmetry of the lattice in the resulting triangular grid, the subdividing diagonal 
is chosen to be along the (+1, + Indirection for square elements in the 1st and 
3rd quadrant [i. e. {k = (k xi k y ) | k x k y > 0}], and along the (+1, — Indirection 
for square elements in the 2nd and 4th quadrant [i. e. {k = (k x , k y ) \ k x k y < 0}] 
of the Brillouin zone. The numerical results shown below were obtained on grids 
with 2N = 2 x 40 x 40 and 2N = 2 x 80 x 80 triangular elements. Comparisons 
between the two grid sizes at a few selected parameter values, including T = 0, 
indicate general agreement of the self-energy results at the 1 — 2%-level or better. 
The chemical potential \i is adjusted so that the electron concentration 44 



equals the specified input value. Note that the prefactor 2 arises from summation 
over the electron spin and f(uj) = l/(e /3uj + 1) denotes again the Fermi factor. 
The self-consistent quasi-particle band Er is then obtained from the solution 




k 



(28) 



of 44 




(29) 



Er =0 . 

kp 



(30) 



The mass enhancement factor on the Fermi surface is 



Z(k F ) ^ 1 -d LU T,'(k F ,uj = Qi). 



(31) 



Corrections to Z due to 44 <%£' are ~ 5 — 
henceforth neglected. 



15%, for the parameters studied, and 
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III. RESULTS AND DISCUSSION 



A. Selfenergy and mass enhancement 

Before discussing the numerical results obtained with the above described 
method, we will first try to obtain a qualitative analytical understanding of the 
self-energy and mass enhancement results obtained at low excitation energies and 
temperatures \u\, T <C fi f . 

For T <C O t , the mass enhancement Z is dominated by the low-energy (~fit) 
lattice tunneling excitations. To estimate this, let us first consider the tunneling 
contribution to the potential, V t , in Eq. (15). By isolating the low-T tunneling 
term, that is, the contribution arising from the excitation i = — > k = 1 from the 
ground- to the 1st excited state, we get 

F/V) = -n\(i; 1 \u\^o)\ 2 [5(nt-w)-S(nt + w)] , (32) 

using fit = t]i — r]o and the fact that the Boltzmann weight filters out the ground- 
state in Eq. (15) at T < fi t , i. e., z^ 1 exp(-/3r7 t ) = 5 ij0 . 

To estimate the transition matrix element in Eq. (32) we introduce the "single- 
well" basis states 

1>±{u) = -±= [Va(«)±V>2(«)J • (33) 

Assume that the overall phase of both ipo{u) and ipi{u) has been chosen such 
that they are both real- valued and both positive for u ~ u + , that is both ifj L (u) 
have positive amplitude in the "right" well of their double well potential w(u). 
Clearly then ifj+(u) must be localized with most of its probability weight in the 
"right" well of w(u), near u = u+, and likewise, by parity, ip-(u) = t(j + (—u) in 
the "left" well, near u = U-. At sufficiently high tunneling barrier, i. e. roughly 
when fi t <^ Q' h = i]2 — i]oi the spatial overlap between •0+ and ifj- in either well 
becomes negligibly small with \ip + (u)\ 2 and \ip-(u)\ 2 sharply peaked near u + and, 
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respectively, near u_. We can thus estimate their -u-matrix elements by a simple 
tight-binding-like approximation, namely 

tyMVv) =U a Sa,a> (34) 

for cr, a' G {+, — }, hence 

(V'lMV'o) = = d*/2 . (35) 

Inserting the foregoing (^1 V'o) into Eq. (32) and the resulting V t into 
Eq. (25), we can now calculate the tunneling contribution to E". This is greatly 
simplified by the fact that V t is ^-independent, i. e., the ^'-integral can be rewritten 
as a frequency integral over the interacting electron density of states 

^) =^ E ^-)G"(^) . (36) 

k 

In the limit T^O and \u\ ~ fi t <C Q' h , we get 

Ej'M = -\r t Q (SU - \oj\) (37) 

where 

r t = \p(Q)C 2 <P t . (38) 

The gap-like depletion of spectral weight in Eq. (37), i. e. the fact that 
E"(u;) = for electronic single-particle energies |u;| up to Q t simply reflects the 
threshold behavior associated with the density of states of the tunneling excita- 
tions: At T = 0, an electron or hole quasi-particle with an excitation energy \u>\, 
measured from the Fermi level, cannot loose more than \u\ in an inelastic emis- 
sion or scattering process. Hence the electron or hole cannot undergo any such 
inelastic processes if \oj\ falls below fit, which is the lowest bosonic excitation en- 
ergy available for emission into the lattice system. Also, at T = 0, there are no 
thermal lattice excitations present to be absorbed by the quasi-particle. Hence 
the damping |E"| vanishes for \ui\ < fi t . 



Inserting Eq. (37) into the Kramers- Kronig relation (27), we find T,' t and, from 
that, the tunneling contribution to the mass enhancement 

Z t = -0 W EJ(O) = ^ • (39) 

Since T t is comparable to electronic energies, whereas fit is small compared even 
to typical phonon energies Q' h , we get Z t 3> 1. In fact, in our model, Eq. (6), we 
have very roughly 

fit ~ n h exp(-c t A B /Q h ) (40) 

with c t of order unity. From microscopic calculations of W, we also know that 
A B ~E P in the EP strong-coupling limit. 6 ' 7 Combining these results, we thus re- 
cover quite naturally in our diagrammatic approach the central result from strong- 
coupling-limit polaron theory, namely 6 ' 7 ' 33 that the electron mass enhancement 
increases exponentially with E P /Q h: that is 

Zt^^-expictEp/Qh) . (41) 

The foregoing results, Eq. (40), (41), apply primarily in the strongly anharmonic 
limit Ab 3> fi/j. However, even in the intermediate regime A B ~ fi/i, fit and 
hence Z t will still remain sensitively dependent on A B , and, via fi^ [Eqs (9), (10)], 
on M and d t . 

In the present diagrammatic model, the T — mass enhancement is thus 
caused directly by the narrow, gap-like depletion in the self-energy spectrum dis- 
cussed above. The detailed spectral shape of E" is of course model-dependent. 
However, the gap-like spectral weight depletion, over a very narrow energy range 
~fi t around the Fermi level, and the overall electronic magnitude of |E"| ~r t ^>fi t 
are generic features of the lattice-mediated interaction Vap in the presence of po- 
laronic anharmonicity. Via the Kramers- Kronig relation, Eq. (27), these features 
of E" are directly translated into the large Z seen in E'. 

Another important feature of the tunneling mass enhancement Z t is its strong 
dependence on the isotope mass M of the lattice degrees of freedom. To see this, 



note first that, all lattice- related model parameters such as A B , d t and C and 
all bandstructure and spin fluctuation-related model parameters, such as t, g s , 
kaf, an d must be regarded as being of purely electronic origin. That is, 
in the general Born-Oppenheimer framework, all these parameters do not depend 
sensitively on M. As a consequence, the damping strength r t , from Eq. (38), is 
also M-independent. However, oc M -1 / 2 , as given by Eq. (9), does depend 
on M. In contrast to this rather modest power-law M dependence of f^, the 
tunneling splitting Q, t is rather sensitive to changes in M. The "exponential" 
M-dependence of f^, reflected roughly in Eq. (40), holds only in the strongly 
anharmonic limit where > and fit *C fi^,. However, Q t will retain a 
substantial M-dependence even at weaker anharmonicity where Ab ~ fi/i ~ fit- 
Hence, a strong isotopic mass dependence of Z t arises in our model, via Eq. (39), 
from the fact that the damping strength T t is independent of M, whereas fl t is 
sensitively dependent on M. 

It is of interest to campare the foregoing results from the anharmonic tunnel- 
ing model to conventional harmonic phonon exchange. To make this comparison 
meaningful, we specifically ask what happens if we replace the anharmonic double- 
well w(u), Eq. (6), by a "corresponding" single- well potential 

w h (u) = l -K h u 2 (42) 

whose the harmonic restoring force constant Kh is chosen to coincide with the 
curvature of the double- well minima in w(u), i. e., Kh = 40As/<i 2 as given by 
Eq. (10) and the corresponding harmonic phonon frequency = ^Kh/M) 1 / 2 by 
Eq. (9). We also assume the same EP coupling constant C. 

The derivation of the self-energy contribution and mass enhancement con- 
tribution from this harmonic Einstein phonon exchange model is standard. 30 ' 31 ' 44 
The self-energy spectrum E'^ has the same analytical form as E" in Eq. (37), with 
Q t replaced by 0,^ and T t replaced by a I\, given analogous to Eq. (38) by 

T h = \p{Q)C 2 dl = p{Q)^l h (43) 
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Here the anharmonic tunneling fluctuation amplitude d t in Eq. (38) has been 
replaced by dh = (2Qh/hKh) 1 ^ 2 which represents, in physical terms, the har- 
monic oscillator zero-point fluctuation amplitude. An important point to note 
here is that, again, dt in the double- well model is entirely a property of the Born- 
Oppenheimer potential function w(u), and hence purely electronic in origin, i. e. 
it does not depend on the isotopic mass M. By contrast, dh, and hence Th do 
depend on M. 

Analogous to Eq. (39), the harmonic phonon exchange contribution to the 
self-energy is estimated by 

Z h = -d^' h (0) = T h /n h = p(0)^ = X h (44) 

which is of course just the conventional dimensionless Eliashberg EP coupling 
parameter A/j, written down here for the simple case of a harmonic Einstein model. 
Comparing Eq. (44) to Eq. (41), we thus find that replacing the harmonic Einstein 
oscillator well Wh(u) by the double well w(u) increases the lattice contribution to 
the electron mass enhancement by a factor 

Z h ~ d\ n t l4DJ 

which will be large compared to unity, if the double well is wide, dt dh, and/or 
the tunneling excitation energy is small, ( <C O/,. 

An essential difference between the harmonic phonon and the anharmonic 
tunneling exchange model is their respective isotopic mass dependence which, in 
Z t , is strong whereas, in Zh, is absent, since cancels out in Eq. (44). The pre- 
dicted isotope dependence of Z t can therefore be used as a crucial experimental 
test to establish or refute the participation of large-amplitude lattice tunneling 
excitations in the low-energy electronic properties of the cuprates. Angle-resolved 
photoemission spectroscopy (ARPES) measurements of the electron quasi-particle 
dispersion, the electronic specific heat, the Drude weight in the low- frequency opti- 
cal conductivity and, in the superconducting state, the low-temperature (T <C T c ) 
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London penetration depth are possible candidates for exploring the isotopic mass 
dependence of the electron mass enhancement in the cuprates. Some experimental 
data for the isotopic mass dependence of the penetration depth in one cuprate 
material already exist 47 and appear to indeed show a strong isotope dependence 
of Z. However, we caution that these data 47 are rather limited, have not been re- 
produced independently at the present time, and their interpretation is presently 
still controversial. 48 

Indirect evidence for a strong isotope effect in Z can be inferred from the 
extensively studied isotopic mass dependence of the superconducting transition 
temperature T c in the cuprates, 49 ' 50 quantified in terms of the isotope exponent 
a = — <91og(T c )/<91og(M) In several non-optimally doped cuprates systems, values 
of a in the range 0.2 — 1.0 are routinely observed, even in samples having T c -values 
as large as 40 - 60K. 49 ' 50 

In theoretical studies of purely electronic AF spin fluctuation exchange pairing 
models, giving rise to a d x 2_ y 2 -pairing instability, one finds 51 that the inclusion of 
harmonic phonon exchange suppresses T c primarily via the phononic contribution 
Zh to the electron mass enhancement. As a consequence, the phononic T c sup- 
pression is largely independent of the isotope mass, due to the absence of isotopic 
mass dependence in Zh, and the isotope exponent is small, typically \a\ < 0.05, 
in harmonic phonon models. Large \a\ values, \a\ ~ 0.5 — 1.0, could be produced 
by harmonic phonon exchange only if the EP coupling in the model was made so 
strong that T c is suppressed to unrealistically low values (< 10K). By contrast, 
in spin fluctuation exchange <i-wave pairing models with anharmonic tunneling 
exchange, of the type described in the present paper, large values of a, of order 
1.0 and larger, could be easily reproduced without substantial suppression of T c . 51 
Taken in conjunction with the experimental isotope data, ' 50 the model calcula- 
tions suggests, that the isotopic mass dependence of lattice contribution to Z is 
indeed much stronger than coupling to a harmonic phonon system would allow to 
produce. 



The foregoing conclusions 51 apply to any electronic "excitation exchange" 
pairing model, with added coupling to harmonic phonon or, respectively, anhar- 
monic lattice tunneling degrees of freedom, and they go well beyond the specific 
framework of the AF spin fluctuation exchange mechanism. We caution however 
that the model isotope calculations 51 are fundamentally based on the assumption 
that the pairing is primarily of electronic origin and that it can be described within 
this general perturbative framework of an electronic "boson exchange" theory. The 
latter assumption is by no means firmly established and or accepted at the present 
time. More detailed experimental studies of the possible isotopic mass dependence 
in the low-energy electronic excitation spectrum of the cuprates would therefore 
be exceedingly helpful in establishing or refuting the participation of anharmonic 
tunneling fluctuations in these materials. 

In interpreting such isotope experiments, especially in the normal state of the 
cuprates, it will be important to recognize that the predicted strong isotopic mass 
dependence of Z is limited to the low-temperature regime. This is illustrated 
by Fig. 1 which shows numerical results for the T-dependence of Z, obtained 
from Eq. (31), and also, separately, the anharmonic tunneling contribution to Z, 
denoted by Zap, both calculated at the Fermi wavevector kp along the (In- 
direction, where 

Z AP {k F ) = -0 w £' AP (fc F ,O) (46) 

and Sap is the contribution to E due to Vap in Eqs. (22), (25). The calculation 
was done with a set of model parameters which is representative for the cuprates, 
namely (rij) = 0.75, that is a hole doping concentration of x = 1 — (rij) = 25%, 
t = 0.35eV, 2 , and spin fluctuation parameters 24 g s = 1.06eV, fi s = 0.105eV and 
1/kaf = 3.33a with a lattice constant a = 3.8 A. Also, 6 C = 3.27eV/A, A B = 
75.29meV, dt = 0.257A, and M = 16u for the lattice atomic mass, so that 6 ' 11 
Qt = 10meV, ft' h = 50meV and 0/, = 97.3meV. 

For T <C Q t , we find numerically that Vap, ^ap and Zap are indeed dom- 
inated by the low-energy tunneling excitation i = — > k = 1 [Eq. (15)] in Vap- 
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Hence, Vap, ^ap and Zap are essentially the same as, respectively, the V t , E t and 
Z t discussed above. Our analytical estimate Z t = T t /CL t , with T t from Eq. (38) 
and the density of states p{u> = 0, T = 0) = 0.175t _1 from the inset of Fig. 2, 
gives Z t = 17.6. This must be regarded as reasonable order-of-magnitude agree- 
ment with the numerical T = result Zap = 8.75 in Fig. 1. The analytical 
overestimate, by about a factor of 2, is due to the fact that our estimate for the 
tunneling matrix element |(0i|'u|V ; o)| — dt/2, entering into F t , is an overestimate. 
Eq (34) is strictly valid only in the large-barrier limit > which is not well 
realized here. Another factor affecting the quantitative accuracy of our analytical 
estimates for E" and Z t is the assumption of an approximately constant p(co), on 
cu-scales of order fl t . In the vicinity of the van Hove singularity [see Fig. 2 inset 
and discussion below], this becomes a somewhat crude approximation. 

Along other directions of the Fermi wavevector, we find numerically quite 
similar results. The overall dependence of Z and Zap on the direction of kp is 
weak. For example, for hp along the (1, 0)-direction, Z is about 10% — 20% larger 
than in the (1,1) direction, with quite similar T-dependence. 

The important point to notice in Fig. 1 is that, while Zap accounts for about 
60% of the total Z at low T, Zap is rapidly suppressed when T becomes compa- 
rable to the tunneling energy scale VL t . This happens because, upon raising T, the 
low-energy spectral gap in S^ p , Eq. (37), gets "filled in," since thermally excited 
quasi-particles can now emit and absorb thermal tunneling excitations and thus 
cause damping, that is |E"| > even at oj = 0. 30 ' 31 ' 44 Already for T ~ Q t , most 
of the 20 t -gap in E" and, along with it, most of the polaronic mass enhancement 
Zap has disappeared. 

Above this characteristic "mass suppression" temperature scale, T~O t , Vap 
contributes only a T- and cu-independent damping to E" which is physically quite 
similar to disorder scattering. The mass enhancement Z is then dominated by 
the spin fluctuation contribution Vsf- In fact, as indicated by the dashed line in 
Fig. 1, for T ~ 2Q t , the mass enhancement in the full "anharmonic phonon + spin 
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fluctuation" (APSF) model becomes indistinguishable from that of the pure spin 
fluctuation (SF) model (without Vap)- This T-independent damping contribution 
can also be estimated analytically for T ^> Q t . The result is roughly 



for all u) <C t. Note that this is the same value as found for low T at w above the 
excitation threshold, u> > Q t in Eq. (37). Note also that this damping contribution 
is independent of the isotope mass M, as discussed above. 

The foregoing high-T results in the anharmonic phonon model should also 
be contrasted with the conventional harmonic phonon exchange model. In the 
latter, the mass enhancment contribution Zh becomes thermally suppressed when 
T ~ Vth- However, the T- dependence of the electron quasi-particle damping con- 
tribution |S^| is quite different in the high-T region, since 30 ' 31 ' 44 ' 52 |£^| ~ \ h T 
rises linearly with T. This difference in the T-dependence is generic to the two 
types lattice dynamics models. These different T-dependences arise because in the 
harmonic model at T > Q h , the amplitude square of the fluctuating EP poten- 
tial seen by the electrons, C 2 {u 2 ), increases linearly with T, by the equipartition 
theorem. By contrast, in the anharmonic tunneling model, this amplitude square 
has a much weaker T-dependence, since C 2 (u 2 ) ~ C 2 d 2 /4 is determined, at both 
low and high T, primarily by fluctuations between two wells, not by fluctuations 
within a single well. 

For comparion, we also quote the approximate, Fermi surface averaged an- 
alytical results for the spin fluctuation contributions to the self-energy, T<sf, 
and to the mass enhancement, Zsf- Using standard Fermi surface averaging 
procedures, 30 ' 31 ' 44 we obtain from Eqs. (17)-(19) and (25) for temperatures T <C 
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where 

T s = ^ p (0)gi(F) FS (50) 

and the Fermi surface average of F(k — k' — Q*) from Eq. (19), denoted by (F)fs, 
is of order unity. The low-energy spectral weight depletion in T,g F (u>) extends to 
fi s . At excitation energies uj on the order of the tunneling energy scale fit, with 
fit <C fi s , the spin fluctuation contribution to the damping is 

l^ F (fi*)| = ^r s (fit/fi s ) 2 = ff^r t «r f . (51) 

This is small compared to the lattice tunneling contribution |E^ p (fi t + + )| = 
(7r/2)r f , even though the respective mass enhancement contributions Zsf and 
Zap may be of comparable magnitude. 

At very high temperatures, comparable to the spin fluctuation energy scale 
fi s , Zsf will be suppressed by the thermal excitation mechanism in the spin 
fluctuation system, analogous to the thermal suppression of Zap described above. 
Because of the much higher energy scale fi s , the on-set of the thermal suppression 
of Zsf is much more gradual, but it can already be clearly seen in the (compared to 
fi s ) relatively low temperature regime T ~ 0.2t shown in Fig. 1. Given their large 
fi s scale, this Zsf suppression effect may not be very important in the cuprates, 
up to, say, 1000K. 



B. Coherent vs. incoherent polaron dynamics 

In spite of the large low-T mass enhancement Zap, the anharmonic lattice 
contribution to the damping in the high-T regime is quite modest, of order |E^ P | = 
(7r/2)rt ~ 0.28eV, using the parameter set from Fig. 1 and Eq. (38). Assuming a 
high-T value Z ~ 4, say [see Fig. 1], this translates into an anharmonic tunneling 
contribution to the quasi-particle decay and transport relaxation rates of order 



This is roughly the right order of magnitude magnitude, compared to typical single- 
particle and Drude peak widths measured in the normal state of the cuprates in 
the 200 — 400K range. 12-19 ' 53 Recall here that the quasi-particle decay rates and 
transport relaxation rates, as reflected in the single-particle and Drude spectral 
(HWHM) peak widths, are being renormalized by the 1/Z-factor, along with the 
quasi-particle dispersion E%. 

It is frequently argued that polarons cannot exhibit coherent transport at 
temperatures exceeding some very low "de-coherence" energy scale, set roughly 
by the single polaron bandwidth, of order 8t/Zp in our 2D model, as estimated in 
the single-polaron strong-coupling theory. 33 ' 34 Here, the strong-coupling polaron 
band width renormalization Zp can be easily as large as 10 — 10 4 (or larger) and 
is exponentially dependent on the EP coupling strength and atomic mass, roughly 
as given in Eq. (41). According to this single-polaron, strong- coupling picture, 
polarons are inevitably condemned to incoherent propagation, by hopping-type 
transport mechanisms, 33 when T becomes comparable to 8t/Z P . Our foregoing 
diagrammatic results suggest otherwise: 

The above estimated high-T decay time, from Eq. (52), and the resulting 
mean free path, are still about an order of magnitude away from the Mott-Ioffe- 
Regel limit, i. e. well on the coherent side of that limit. The quasi-particles, 
albeit damped, remain coherent at high-T, due to the fact that the polaronic mass 
enhancement Zap is being suppressed at T ~ Q> t and due to the fact that the 
high-temperature "remnants" of the low-T anharmonic tunneling fluctuations do 
not scatter the quasi-particles very strongly. At least for the specific parameter set 
used above, the scattering from the "remnant" anharmonic fluctuations is by no 
means sufficient to drive the electron quasi-particles into the incoherent / hopping 
type regime. The physical picture that emerges here is that, at T ~ only 
the coherence of the lattice tunneling groundstate - but not the coherence of the 
electronic quasi-particles - is being lost when T reaches the lattice tunneling energy 
scale. 
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How can these two seemingly quite contradictory physical pictures be recon- 
ciled ? The answer, we believe, lies in the parameters. The above-described strong- 
coupling small polaron picture is based on the EP strong-coupling limit where, in 
the underlying microscopic EP model, the EP coupling strength Ep = C 2 / K ex- 
ceeds the bare bandwidth 8t and, at the same time, the tunneling barrier heights 
As of the polaronic Born-Oppenheimer lattice potential W and the polaron bind- 
ing energies (relative to the delocalized carrier state) are large, namely, of order 



In this limit the polaronic lattice fluctuation amplitude dt is of the order 6 ' 7 
d t ~ C/K, hence Cd t ~ Ep, and in our local double- well lattice model, Eq. (6), 
the condition Ep 3> t, via Eq (38), translates into 



We are assuming here p(0) ~ l/8t so that p(0)Ep ~ 0(1). Under these conditions, 
our diagrammatic approach breaks down at T ~ but, nevertheless, the dia- 
grammatics signals its own inadequacy, as the quasi-particle damping, estimated 
from Eq. (47), will indeed approach or exceed the Mott-Ioffe-Regel limit. 

This strong-coupling picture is inevitably correct for low-density polaron sys- 
tems with a nearly empty conduction band and with short-range EP coupling in 
spatial dimensions D > 2. In this type of EP system, the conditions Ap ~ Ep 3> t 
must be satisfied in the polaronic regime because, unless Ep > Ep rit ^ ^> t, po- 
larons will simply not be formed. 35 Here, Ep Tlt ^ denotes the minimum EP coupling 
strength needed to stabilize the polaronic self- localized carrier state against der- 
ealization. In the nearly empty-band situation, Ep rit ^ is of order of the bandwidth 
or larger. 6 ' 35 

However, in the case of the Holstein-Hubbard (and related) models, as applied 
to the cuprates, polaron formation arises in the near 1/2-filled band limit and the 
polaronic carriers are dopant induced holes in a strongly correlated Mott-Hubbard 
insulator electron background. Extensive microscopic calculations 6-9 have shown 



Ep. 




(53) 
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that the presence of AF spin correlations in the near 1/2-filled Hubbard electron 
system drastically reduces the threshold for polaron formation in such a system, 
roughly to Ep Tlt ^ ~ 2 — 3t in realistic parameter regimes. Furthermore, for Ep 
near Ep Tlt \ the lattice potential barriers separating different polaronic minima 
can become quite small, so that A B <C t or even A B ~ &h and also Cd t < Ep. In 
this limit, we get 

r t < \p{G)E 2 P ~ t (54) 

In this " weakly bound polaron" regime, a diagrammatic treatment based on a 
coherent quasi-particle picture is likely to become applicable, even for T ~ tt t . On 
the other hand, the strong-coupling expansion, based on the assumption Ep ^> t, 
is likely to fail here. 

It is conceivable that microscopic EP models may exhibit a cross-over from 
the strongly bound regime, Eq. (53), to the weakly bound regime as a function of 
increasing doping and/or increasing temperature. Such a cross-over can be driven 
by mutual screening of the attractive polaronic EP potential wells. This screening 
occurs at finite polaron density when nearby polaronic self- localized wavefunctions 
and EP potential wells begin to overlap. It has the effect of lowering the height 
A B and the width dt of the polaronic inter-site tunneling barriers in the lattice 
potential W. This barrier screening effect has recently been demonstrated by 
microscopic calculations for the Holstein- Hubbard model. 7 Thermal de-trapping 
of the polaronic carrier from its EP potential well will also tend to enhance this 
barrier screening effect by incrasing the population of unbound carriers. Hence, 
raising the temperature can have qualitatively the same effect as raising the doping 
level, namely, to push the system towards weaker polaronic binding energies and 
lower barrier heights. It will be interesting to investigate whether such a cross- 
over is realized in at least some of the cuprates which, as a function of doping and 
temperature seem to show evidence of hopping-type conductivity at low doping 
and low T, but exhibit apparently coherent transport at larger doping and/or 
higher T. 
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The partially ionic character of the cuprates and resulting spatially extended 
Frohlich-type EP couplings are another factor which may contribute to the exis- 
tence of a "weakly bound" polaron regime. Such long-range EP couplings have the 
effect of lowering the overall threshold Ep rit ' by allowing "large" (i. e. spatially 
extended) polaron formation to occur even in dimension D > 2. 35 ' 54 Extensive 
phenomenological studies have so far considered only large bi-polarons and only 
in the context of continuum models. 35 ' 54 It will be of considerable interest to in- 
vestigate such longer-range EP coupling with discrete crystal lattice models and 
to explore the cross-over from large / weakly bound to strong-coupling / strongly 
bound polaron behavior. 



C. Density of states and single-electron spectra 

In Fig. 2, we show the resulting density of states, p(u>), calculated for the 
same parameter set as in Fig. 1, from Eq. (36), at T = and at T = 0.029t = 
The most prominent feature in Fig. 2 is the narrow, strongly T-dependent peak 
in close proximity to the Fermi level uj = in p(u), as shown by the full line in 
Fig. 2. 55 At T = 0, we can show that p{oS) near the Fermi level is essentially just 
a narrowed image of the non-interacting density of states 44 p^°\ui) for the same 
band filling, that is, for |u;| 

p{u)^p^\Zu) (55) 

— » 

with Z denoting the /c^-average of Z. The "resonance" peak in Fig. 2 is thus 
primarily a narrowed image of the non-interacting 2D van Hove singularity and 
its peak position uj v h s relative to the Fermi level u = is roughly given by 

(o) 

UvHs — 2 

where oj^Hs denotes the position of the van Hove singularity, relative to the Fermi 
level, in the corresponding non-interacting system at the same band filling. 
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A rough derivation of this result, Eq. (55), can be sketched as follows: First, 

— * _ — * 

neglect the k dependence of E around the Fermi surface and replace E(/c, u) by an 
appropriate Fermi surface average E(u;). 30 ' 31 ' 44 Second, neglect the small imagi- 
nary part E" and approximate its real part by the standard "Fermi liquid" ap- 
proximation 

E'H = (1 - Z)oj + O(uo 2 ) (57) 

where 2 = 1 — <9 W E(0). Third, insert this into Eqs. (26) and the resulting G" into 
Eq. (36) to arrive at Eq. (55). 

Note, in Fig. 2, that the logarithmic divergence at the van Hove singularity 
has been slightly smeared out, even at T = 0, by the damping (|E"| >0) at finite 

— * 

u and by the k dependence of E along different directions on the Fermi surface. 

— > 

These damping and /c-dependence effects are neglected in the derivation of Eq. (55). 
Also, as shown in the inset of Fig. 2, the resonance peak is rapidly suppressed, 
along with Z(T), when T becomes comparable to Q t . This is again caused by the 
strong thermal damping at lv = and by the thermal suppression of Z(T), when T 
bcomes of order Q t or larger. Because of the thermal damping effects, the density 
of states peak in the "APSF" model (including Vap) is suppressed faster with 
increasing temperature than in the pure "SF" model (without Vap-, see dashed 
line in Fig. 2 and inset); the APSF model generates an additional strong damping 
already when T ~ Q t which is absent in the SF model. 

A careful study of the doping dependence of the density of states peak near 
lo = reveals that, while substantially contributing to this peak, the van Hove 
singularity in the non-interacting density of states is not necessarily required 
to generate such a p(w)-peak. Even for a completely flat p(°\u), e. g. at band 
fillings where the van Hove singularity is far away from Fermi level lo = 0, it is 
possible to obtain a pronounced peak structure in p(u>) near the Fermi level. The 
general conditions under which this happens are that the self-energy spectrum 
\Y,"(k F , lo)\ on the Ferrmi surface kp must (i) be very small (or zero) at to = and 
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(ii) rise very rapidly to large values on some very low frequncy scale, denoted by 
fi, such that 



Note that these conditions can be satisfied in our APSF model at temperatures 
T C ( , where Q t represents the low-energy scale and it is therefore possible 
to get X"(Ot + + ) ~ T t 3> By contrast, these conditions can, realistically, 
not be satisfied in harmonic phonon exchange models, as exemplified by the har- 
monic Einstein model discussed above. In the latter, we have, from Eqs. (43) 
and (44), £^(0^ + + ) = (ir /2)\ h Q h: which is comparable to the phonon energy 
scale for typical Eliashberg A^-values of order unity. Hence, EP-induced density 
of states peak effects are another important distinguishing feature between har- 
monic phonon and anharmonic tunneling exchange. We will now describe how this 
density of states peak effect arises. 56 

Under the conditions of Eq. (58), the quasi-particle peak in G" is rapidly 
damped out with increasing quasi-particle energy which, in essence, has the effect 
of pushing spectral weight away from the Fermi level for finite u, of order O. This 
spectral weight removal, due to strong damping at u> ~ O, occurs in addition to 
the spectral weight removal caused by the mass enhancement factor Z. In contrast 
to the latter, the former spectral weight removal is not being compensated for by 
the reduced quasi-particle dispersion, when one averages G" over all k to calculate 
p(ui). At \oj\ <C O, on the other hand, spectral weight in G" is removed only by 
the Z-effect which, in p(u), is being compensated for by the reduced quasi-particle 
dispersion. As a consequence, spectral weight will be removed from p{ui) for finite 
uj, of order Q, but not at very low uj, well below O. It is important to realize that, 
therefore, the resulting peak in p(u) at u = does not arise from an enhancement 
of p(oj) at the Fermi level, but rather from a suppression of p(oj) off the Fermi 

— * 

level. In fact, neglecting the /c-dependence of S, the density of states right at the 
Fermi level is essentially unchanged by the interaction, that is, from Eq. (57) 



z"(k F ,n)\ ~t > o, |e"(£ f ,o)| . 



(58) 



P(0) = P (0) (0). 



(59) 



Notice also that a p(u;)-peak, generated solely by the foregoing damping- 
induced spectral weight removal mechanism from a flat non-interacting density of 
states (i. e. from a p(°)(u;) = const at \u>\ ~ fi), will be centered right at the Fermi 
level lo = 0. As shown in the inset in Fig. 2, this is not the case in our T = model 
calculation. Here, the narrowed image of the van Hove singularity, centered at an 
^vHs — 0.02t, is below Q t and superimposed on the damping-induced spectral 
weight suppression which occurs mainly for \u\ > fit in our model, due to the 
almost completely gapped E" at \u\ < fit, Eq. (37). However, if we change the 
band-filling so as to move the van Hove singularity further away from the Fermi 
level, then the narrowed van Hove image moves out to larger uj v h s and gradually 
disappears when u v h s becomes comparable to fit while, at the same time, the 
damping-induced peak at u = emerges, abeit generally with less pronounced 
peak height. 

We emphasize the foregoing especially in relation to photoemission studies 
of the /c-integrated density of states and its doping evolution in the cuprates. 
The foregoing considerations suggest that one needs to be very careful in the 

— * 

interpretation of observed /c-integrated density of states peaks. Such peaks need 
not necessarily be a manifestation of van Hove singularities, but may rather be 
caused by other effects, such as strongly frequency dependent damping near the 
Fermi level. 

Also shown in Fig. 2, by the dot-dashed line, is the result of a non-self- 
consistent 1st order SF model calculation for p(u). Here, the fully dressed internal 
single-electron Green's function in Eqs. (22), (25) is replaced by the bare Green's 
function at the same band filling, that is, written in the Matsubara domain, with 

G<°>(£,i„) = ^_. (60) 
iv — e T 

k 

Here ej-p denotes the non-interacting band, obtained from Eq. (24) with the non- 
interacting value of the chemical potential fj,(°> at the same band filling (rij)^ = 
0.75. 
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Similar 1st order calculations were first reported for this SF model in Ref. [24]. 
An important prediction from these early 1st order model calculations was the 
demonstration that the spin fluctuation exchange causes a pronounced "pseudo- 
gap" to open up in the density of states near the Fermi level. This pseudo-gap, 
i. e. a strong suppression in the density of states near to = 0, is indeed also 
present in our 1st order calculation, with the spectral weight depletion extending 
approximately from u> ~ up to u ~ 3.8t in Fig. 2. However, in our self-consistent 
SF calculation (using a fully dressed G as the internal electron line), the pseudo 
gap has largely disappeared, as shown by the dashed line in Fig. 2. In the w-region 
of the gap, where the 1st order p{uj) is being depleted, the self-consistent p{uj) is 
almost flat. Hence, it is not quite clear whether the original claim of pseudo-gap 
generation in AF spin fluctuation exchange exchange models is really tenable. At 
the very least, this spin fluctuation induced pseudo-gap effect, if present at finite 
doping, appears to be quite model dependent, that is, very sensitive to the level of 
approximation employed in the calculation and, possibly, to the choice of model 
parameters. 

In Fig. 3, we show the single-hole excitation spectra 44 

P(k 1 u;)^n- 1 \G"(k^)\f(u;) 1 (61) 

with the Fermi factor f(u) = l/(e Pu + 1), for T = 0.029t = O t and several Appoints 
along the (1, 1) direction. At all k, P(k,u) exhibts a strong continuous "incoher- 
ent" background which is almost constant, down to u> ~ — 2 t, and extends up to 
the Fermi level w = for /c = As k — > kp(l, 1), the expected quasi-particle 
peaks at energies E% < emerge from this background, sharpen, and disperse 
towards E% = 0. These results are quite similar to typical angle-resolved photoe- 
mission spectroscopy (ARPES) data obtained in the cuprates at T > T c . 12-19 ' 57 
Notice that, at the temperature T = 0.029t used here (i. e. T ~ 110K, assuming 
t= 0.35eV), the lattice contribution to the mass enhancment Zap is already no- 
ticeably reduced by the thermal Z-suppression. Hence, except for the additonal 
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line broadening caused by £^ P , the shapes and the dispersion of the quasi-particle 
peaks are already very similar to those in the purely electronic SF model, with- 
out coupling to lattice degrees of freedom. We emphasize that the normal-state 
ARPES line shapes and quasi-particle dispersion in the cuprates are measured at 
typical T>T C ~ 80 — 120 °K. Even though coupling to lattice tunneling excitations 
may be important in shaping the quasi-particle dynamics of the cuprates at low 
temperatures, the tunneling excitations do not necessarily have any pronounced 
signature at normal state temperatures which may well be in the regime T ~ Q, t . 

The quasi-particle energy band -Eg, shown in Fig. 4, also resembles typical 
cuprate ARPES data near k F . 19 In particular, the very flat ("heavy mass") por- 
tions in Er near kp along the (1, 0)-direction are very similar to recently reported 
quasi-particle dispersion curves along certain /c-directions in the cuprates. 19 In our 
model, these very "heavy" sections in E% are due to the close proximity of /c P (l, 0) 
to the non-interacting band saddle point, the van Hove point (tt, 0), with substan- 
tial additional band "flattening" caused by the mass enhancement Z. Note here 

— * 

that Eq. (29) has multiple solutions for E% near kp, as shown in Fig. 4. The E%- 
branch crossing the Fermi level E% = corresponds to the quasi-particle peak in 
G"(k,u). Other branches do not necessarily produce well-defined separate peaks 
in G" ', because of damping. 



IV. CONCLUSION 

In summary, we have proposed a highly simplified anharmonic lattice potential 
model to study the effect of polaronic tunneling excitation and anti-ferromagnetic 
(AF) spin fluctuation exchange on the quasi-particle dynamics in a 2D single band 
model. The effects of both polaronic lattice tunneling and AF spin fluctuation ex- 
change are treated in a self- consistent Migdal-type single-exchange approximation 
for the single-particle self-energy. At temperatures T and excitation energies u> 
which are well below the anharmonic tunneling excitation energy scale fit, the 



single-electron self-energy is dominated by the exchange of lattice tunneling exci- 
tations which give rise to a large and strongly isotopic mass dependent electron 
mass renormalization Z. 

With increasing T, the lattice tunneling contribution to the electron mass 
enhancement, along with the isotope effect of Z, are rapidly suppressed and they 
both essentially disappear when T exceeds In the high temperature regime 
T ~ 2Q t , the electron mass enhancement Z is completely dominated by the AF 
spin fluctuation exchange; the lattice tunneling excitations contribute only to the 
electron quasi-particle damping |E"| and this damping contribution is largely T- 
independent and independent of the isotope mass. This high-temperature damp- 
ing effect is therefore quite similar - and experimentally indistinguishable from - 
conventional disorder scattering. 

Surprisingly, this anharmonic lattice contribution to the quasi-particle damp- 
ing at high temperatures T ~ Q t can be quite small compared to the electronic 
bandwidth, even if the anharmonic lattice tunneling fluctuations produce a very 
large mass enhancement at low T. The relatively small damping at high T allows 
the electronic quasi-particle propagation to remain coherent at high temperatures, 
despite of the polaronic character of the underlying lattice dynamics. We have de- 
scribed in detail the "weakly bound polaron" conditions under which this behavior 
may occur, within the framework of a microscopic Holstein-Hubbard-type model; 
and we have contrasted it with the fundmentally incoherent high-T behavior in 
"strongly bound" polaron systems. 

In comparing the isotopic mass dependences of the electron mass enhance- 
ment Z at low T in the anharmonic tunneling exchange model versus a conven- 
tional phonon exchange model, we find that the two types of models exhibit fun- 
damentally different behavior: The Z contribution in the anharmonic tunneling 
exchange model is isotope dependent, whereas, in the harmonic phonon model, it 
is not. Furthermore, in contrast to the T-independent damping in the anharmonic 
tunneling exchange model at high T, the high-T damping contribution in the har- 



monic phonon model increases linearly with T, In this context the demarkation 
between "low T" and "high T" is given, respectively, by the tunneling energy scale 
Qt in the former and by the harmonic (Einstein or Debye) phonon energy scale in 
the latter model. 

Based on the foregoing results, we suggest to explore the isotopic mass de- 
pendence of the low-T electron mass enhancement experimentally, in order to 
establish whether or not anharmonic lattice tunneling excitations contribute in a 
significant manner to the low-energy electronic properties in the cuprate mate- 
rials. Possible candidate experiments for pursuing such investigations are quasi- 
particle energy dispersion curves, as measured in angle- resolved photoemission 
spectroscopy (ARPES); Drude spectral weights, as observed in the optical con- 
ductivity; electronic specific heat; and, in the superconducting phase, the T — > 
London penetration depth. 

In the /c-integrated density of states, low-T tunneling excitation exchange can 
produce a narrowed image of the near-Fermi-level features of the non-interacting 
density of states, such as, for example, a narrowed van Hove singularity peak. At 
finite temperatures T ~ fi t , the van Hove peak is completely suppressed due to the 
thermal damping. Likewise, the quasi-particle dispersion Er is largely unaffected 
by the anharmonic tunneling excitation exchange and the single paricle spectral 
function, as observed in ARPES experiments, is affected only by the T-independent 
broadening, when T ~ Q t . 

We note that our results concerning the low-T narrowed imaging of the van 
Hove peak and its thermal suppression are generic to any strong-coupling low- 
energy boson exchange model, as demonstrated by our results for the purely elec- 
tronic spin fluctuation exchange model (see Fig. 2). In boson exchange models 
of the type discussed here, 58 the rapid thermal suppression of the narrowed van 
Hove peak in p(u>) is ubiquitous and entirely in line with the finding that van 
Hove singularities alone are not sufficient to bring about or substantially enhance 
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high-T c pairing instabilities. Our results thus strongly support the notion 58 that 
a conjectured "van Hove scenario" 59 of high-T c superconductivity does not exist. 

To the extent that the cuprates can be described in terms of this type of 
low-energy "tunneling boson" exchange picture, we emphasize that the observed 
normal-state mass enhancements 12-19 Z(T) ~ 2 — 4 may substantially under- 
estimate the true bosonic coupling strength 44 , A = Z(T = 0) — 1. As shown 
here, Z(T) <C Z(T = 0) if the normal state measurement temperature T (e. g. 
T>T C ~ 80 — 12QK in the cuprates !) exceeds the relevant boson energy scale O. 
Also, for T<^cO, the full Z(T = 0) will be observable only within a very narrow en- 
ergy range around the Fermi level, which could well be below current limits 
of resolution. 12-19 In high-T c systems where Q t ~ T c for T<^T C , it needs to be 
explored to what extent the full normal-state value of Z(T = 0) is observable in the 
presence of a superconducting gap A with A~0. Finally, we note that a strongly 
T-dependent Z has important implications for the T-dependence of quasi-particle 
lifetimes r, since l/r=|E"|/Z is renormalized by Z. 

Given sufficient energy resolution, the full T = normal-state mass enhance- 
ment should be observable in systems where T c <^fi, at temperatures T well below 
O. From the T-dependence of Z(T), the relevant boson energy scale can then be 
established. Also, the possible observation of the p(c<;)-resonance peak should be 
explored. The resonance is observable by either direct or inverse angle-integrated 
photoemission spectroscopy, depending on the resonance peak position uj v h s rel- 
ative to the Fermi level lo = 0. The resonance can also cause a low-energy gap 
structure, on an energy scale of order |c<;„h s |, in the two-particle electronic re- 
sponse functions. 

To put our model results in a more general perspective, we note that the 
specific shape of the spectral weight depletion around the Fermi level u = in the 
self-energy [see Eq. (37)] is obviously model dependent. However, the qualitative 
feature of its existence is characteristic of the low-energy behavior of any Fermi 
liquid. What is unusual about this Fermi liquid behavior in the lattice tunneling 



self-energy studied here, is the fact that the Fermi liquid character begins to break 
down already at a very a low excitation energy (u) and temperature scale, set 
by the tunneling energy Q t . This can, but need not be lower than all other 
electronic or lattice energy scales, because of the great sensitivity of tunneling 
matrix elements to small changes in the tunneling barriers. Thus, the tunneling 
energy scale (i) could vary substantially between different materials; (ii) within 
the same parent material, it could vary substantially as a function of doping level; 
and (iii) it could well be lower than either the harmonic phonon or even the 
superconducting T c scale in the cuprates. 

These last features of our model are generic to the low-energy physics of 
strongly correlated polaronic electron systems; they do not depend on the model 
simplifications we have introduced here. What is quite likely to change in a more 
realistic polaron model is the structure of the anharmonic lattice dynamics, as re- 
flected e. g. in the ui- and ^-dependence of Vap- This, in turn, will have important 
ramifications for the details of the T- and (^-dependence of the self-energy in such 
systems. 

Lastly, we should also emphasize that the Hubbard-type local Coulomb cor- 
relations play a central role in such a polaronic picture, for two reasons. Firstly, 
the local Coulomb repulsion is important to prevent the formation of bi-polarons 
and hence, to preserve the essentially fermionic character of the quasi-particles. 
Secondly, the Coulomb correlations are the origin of the AF spin correlations in 
the nearly ^-filled regime. As discussed previously, 6-9 the AF spin correlations 
are essential to allow doping induced carriers to form polarons at much weaker 
EP coupling strengths than would be required in weakly correlated systems. As a 
consequence, the resulting polaronic system can have much lower tunneling barri- 
ers and hence, reach much larger O t -scales than one would estimate for weakly / 
un-correlated small-polaron systems. 

Note that the latter is crucial if one wants to apply such a polaronic picture 
to the cuprates. On the one hand, one will need low tunneling energy scales in 
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order to explain the breakdown of Fermi liquid behavior down to very low temper- 
atures, down to, say, ~ 10K, in some cuprate systems; on the other hand, one will 
need large tunneling energy scales to to explain the unusual T-dependences of the 
quasi-particle and transport dynamics which are found over very large tempera- 
ture scales, up to ~ 100 — 1000K. A lattice excitation spectrum covering several 
orders of magnitude of excitation energy will therefore be essential. This, again, 
emphasizes the importance of more realistic models of the polaronic lattice dy- 
namics and will require further study. 
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FIGURE CAPTIONS 



FIG. 1. Mass enhancement Z = 1 — d C0 'E'(kF,uj = 0) vs. temperature T at the 
Fermi surface point &f(1, 1) =0.44(-7r, 7r) in the APSF model, with V = Vap + Vsf 
in Eq. (25), (upper full line) and in the purely electronic SF model, with V = Vsf 
in Eq. (25) (lower full line). Also shown is Zap = —d UJ T, , A p(<jj = 0) vs. T in 
the APSF model (dashed line), where T,' AP is the V/ip-contribution to £' in the 
APSF model. The model parameters are given in the text in Sec. III. A, following 
Eq. (46). 

FIG. 2. Density of states p(uj) vs. energy uj at temperature T = Q t - APSF (full 
line) and SF (dashed line) are the results, for the respective models, as defined in 
Fig. 1. The results labeled SF1 (dot-dash line) are results of non-self-consistent 
1st order calculations for the SF model, as descibed in text. Inset: APSF and SF 
results for p(u) vs. ui near cu = at T = Vt t and T = 0. The other model parameters 
are the same as used in FIG. 1. 

FIG. 3. Hole excitation spectra P(k, u) vs. energy u in the APSF model [with 

V = Vap + Vsf in Eq. (25)] for k = (k x , k x ) at /c x -values indicated, at temperature 
T = Q t . Spectra are off-set by successive vertical shifts, with dashed lines showing 
zero intensity. The other model parameters are the same as used in FIG. 1. 

FIG. 4. Quasi-particle bands vs. wavevector k in the APSF model [with 

V = Vap + Vsf in Eq. (25), full line] and in the non-interacting system [V = 0, 
dashed line], along several directions in the 2D Brillouin zone, both at the same 
electron concentration (nj) =0.75, at temperature T = Vt t . The other model pa- 
rameters are the same as used in FIG. 1. 
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